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Abstract. In this paper we derive approximate quasi-interpolants when the values of a 
function u and of some of its derivatives are prescribed at the points of a uniform grid. As a 
byproduct of these formulas we obtain very simple approximants which provide high order 
approximations for solutions to elliptic differential equations with constant coefficients. 
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1 Introduction 

The method of approximate approximations is mainly directed to the numerical solution 
of partial integro-differential equations. The method provides simple formulas for quasi- 
interpolants, which approximate functions up to a prescribed precision very accurately, 
but in general the approximants do not converge. The lack of convergence, which is not 
perceptible in numerical computations, is offset by a greater flexibility in the choice of ap- 
proximating functions. So it is possible to construct multivariate approximation formulas, 
which are easy to implement and have additionally the property that pseudodifferential op- 
erations can be effectively performed. This allows to create effective numerical algorithms 
for solving boundary value problems for differential and integral equations. 
Approximate quasi-interpolants on a uniform grid are of the form 

Mn{x) = V--/^ J2 ^ i'-r^) ^(^^•) (1) 

with positive parameters, "small" h and "large" T), and the generating function Ti is 
sufficiently smooth and of rapid decay. Their properties have been studied in a series of 
papers [4-7] (cf. also the monograph [11]). It has been shown that the quasi- interpolant 



approximates smooth functions with 



N-l 

\Muix) - u{x)\ < e J2 {hVvf^\d''uix)\ + c{hVvf\\VNu\\L^ , (2) 

as long as 7i is subject to the moment condition 

= 5|„|o , 0<|a|<A^. (3) 

In ([2]) the constant c depends only on 7i and e can be made arbitrarily small if the 
parameter D is sufficiently large, so that one can fix T> such that in numerical computations 
Mu approximates with the order 0{h^). The construction of simple generating functions 
satisfying ^ for arbitrary has been addressed in [9], whereas [3,10] extend the results 
to the case of nonuniform grids. 

In this paper we study more general quasi-interpolation formulas with values of a 
function u and of some of its derivatives prescribed at the points of a uniform grid. More 
precisely, we consider approximants of the form 

Mu{x) = n f ^^"l Q{{-h^)d)u{hj) , (4) 

where Q(t), t E W^, is a polynomial with deg Q < N. 

We establish estimates of the type ([2]) if the generating function 7i and the coefficients 
of Q are connected by suitable conditions (see (jl4p ). It is shown in particular, that 
for an arbitrary polynomial Q there exists TC such that the approximate Hermite quasi- 
interpolant satisfies the estimate ([2]). On the other hand, the same is true if 7i is the 
Gaussian or a related function and Q is chosen suitably. As a byproduct we obtain very 
simple approximants which provide high order approximations to solutions of the Laplace 
equation or other elliptic differential equations. 

Estimate ([2]) is proved in Section [2] for ^ with TC and Q connected by conditions 
(|14p . In Section [3] we describe a simple method for the construction of TC satisfying the 
conditions (jl4p for a given polynomial Q. This result is applied in Section [J] when we 
consider some examples of Hermite quasi-interpolants. In the particular case TC{x) = 
^-n/2|-jg^ ^^-i/2g-<_B x,x> ^ with B = {hj} Symmetric and positive definite real matrix, 
we obtain the following quasi-interpolant of order 0{{h^/V)'^^^) 



where B is the second order partial differential operator Bu = bijdidjU. In Section 

Owe use Mu for the approximation of solutions of the equation Bu = 0. If n satisfies the 
equation in then Mu has the simple form of the quasi-interpolant of order 0{{h^/V)'^) 
but gives an approximation of exponential order plus a small saturation error. The same 
approximation property holds for functions which satisfy the equation Bu = in a domain 
Q C M". If the function u is extended by zero outside O then Mu simplifies to 

Mn(x) = ^^^^^ y n(/im)e-<^"'(--'^-)'--'^™)/('^'^). 
^ ^ (7rP)"/2 ^ ^ ' 

We obtain that, for any e > 0, M u approximates pointwise ti in a subdomain Q' C with 

2M-1 

\Mu{x)-u{x)\ <e Y {hVv)\'^\\df^u{x)\+CB{hVvf^Cu, 



if we choose T) and h appropriately; the constant is independent of n, h^V, M and Cu 
depends on u. 

In Section [6] we show that the Hermite quasi-interpolant (jH) gives the simultaneous 
approximation of the derivatives of u. If d^^Ti exists and is of rapid decay, for any function 
u G with L>N + \(5\, the difference df^Mu{x) - d^uix) can be estimated by 

L-l 

\d^Mu{x) - d^u{x)\ <e + ci ihVv)^\\VN+\p\u\\L^, 

|7|=0 

with ci independent of u, h and V. 



2 Quasi-interpolants with derivatives 

In this section we study the approximation of a smooth function x £ R", by the 

Hermite quasi-interpolation operator 

Mu{x) = n ( ^-^ \ Q (^-hVVd^ u{hm) (5) 

where Q{t), t S M", is a polynomial of degree at most — 1 

N-l 

Q{t) = Y «7*^' * ^ with G ffi and ao = 1 , (6) 

|7|=0 

d = di . . . dn, and W is a sufficiently smooth, rapidly decaying function. Then 

Our aim is to give conditions on the generating function TC such that d?]) is an approx- 
imation formula of order ©((/iV^)^) plus terms, which can be made sufficiently small. 

Suppose that u £ C^(R") n L°°(M"). Taking the Taylor expansion of d'^u at each node 
hm leads to 



d^u{hm)= Y , d''+^u{x) 

\a\=0 

E(hm — x)" / , ^ 

^ Uo,+^{x, hm), < |7| < iV, 



(8) 



with 

1 



f/„,„) = A./,-'a"„(,. + (i-,).).„,N = iv^ (9) 



We write Mit in the form 

N~l 



MU{X) = Y i-flVVy^^d^uix) J2 ^C7a{l,V,n)+nH,N(.^ 



=0 o</3 



with the periodic functions 



and the remainder term 



|/3|=7VO<o</3 ■ meZ" 

Therefore by using the definition of a a, 
Mu{x) - u{x) = u{x) (ao{^,V,n) - 



n 



hVv J V hVv J 



■ 



N~l 



|/3|=1 Q</3 



(10) 



Let us introduce the functions 



£a{L n) = V,n)- j x^nix) dx, O <\a\ <N -l. 



Then we have 



ao{i,V,H)-l = £o{i,V,H) + (^j H{x)dx-l 
and for = 1, . . . , iV - 1 



Theorem 2.1 Suppose that Ti. is dijferentiable up to the order of the smallest integer 
uq > n/2, satisfies the decay condition: there exist K > N + n and Cp > such that 



\d^nix)\<CB{i + \x\)-^, o<|/3|< 



no 



(13) 



and the conditions 



J n{x)dx = 1 

j x'"H{x)dx = Q , |/?| = l,...,iV- 1. 



(14) 



Then for any e > there exists V > such that, for all u G W^(M") n C^(M") the 
approximation error of the quasi-interpolant ([7]) can be pointwise estimated by 



N-l 



\Mu{x) - u{x)\ < e Y {hVvf^ldfuix)] + cihVv)'' \\V nu\\l^ (15) 



with the constant c not depending on h, u and T>. 



Proof. Under conditions ()14p the relations (llip and (1121) simplify to 



«</3 a</3 



Moreover, the assumptions (fT3]l on H. ensure that (see [10]) 



and 



as D ^ oo. 



J2 \d"Tn{Vvu)\ 

;/eZ"\0 

Hence for any e > there exists D such that the following estimates are valid 



^aaiC,V,n)\ < ^ l^\£^iC,V,n)\ <e, = 1, . . . , AT - 1. 
■"^ a! I a! 



|ao(^,P,W)-l| = |fo(e,^?,W)| <e, 



a</3 a</3 

The next step is to estimate the remainder term TZh^^. Since 

1 



\U^{x,y)\=N 



s^-^d'^u{sx+{l - s)y)ds 



we obtain 



< Wu\\l^ 



\P\=N 



where 



0<a<P 



^ — m \ " — m 
ri 



(16) 
(17) 



In view of the decay condition there exist constants such that, for T) sufficiently large 
( [10]) 

||p„(-,P,W)||l^ <c„, 0<|a|<iV. 



Hence we obtain 



IL°° 



(18) 



m=N 



with a constant c not depending on /i, T> and u. By (jlOp . this leads together with ([16 
and p!7|) to the estimate ([15]). 



3 Construction of generating functions for arbitrary N 

Here we describe a simple method for the construction of generating functions TC satisfying 
the conditions for arbitrary given with ao = 1. Let us denote by ^ = {^4^^} the 
triangular matrix with the elements 

J a/3-a a<P, ifli I I n at 
Aaf3=<^ . |/3|,|a| =0, ...,iV-l. 

L) otherwise 



The dimension of A is (A^ + n — 1)!/((A^ — 1)! n!). Since det^ = 1 there exists the inverse 
-1 

1/3 



matrix A ^ = {-^Qfl^^} and (jl4p leads to the following conditions 



x''n{x)dx = alA^jK |a| =0,...,A^-1. (19) 

These conditions can be rewritten as 

9°(^W)(0) = a!(-27ri)l"US,;^\ |q| = 0, . . . , iV - 1. 
Let us assume (see [9]) 

1 d 



where VN{t) is a polynomial of degree less than or equal to — 1 and is a smooth 
function rapidly decaying as |x| — > oo with J-rj[{}) ^ 0. Conditions (jl9p give 

a°(P;v(A)^7?(A))(0) = (-27ri)Ha!4;^\ |a| = 0, . . . , - L 
We choose Vn as the Taylor polynomial of order {N — 1) of the function 

Since 

5^Q(0) = J]4;i)(-27ri)l"l-^5'^-"(^r/)-i(0), 
where we use the notation 

a/3-(^^)-i(o) = a'3-(_L_)(o), 

we obtain 

|/3|=0 a<l3 ^'^ 

Therefore the equations 

a-(Pjv(A)7-7?(A))(0) = a°(Q(A)^7KA))(0) 

= 4-^)(-2vriA)^)(0) = (-2^0l"la!4;^) 

|7|=0 

are valid for alla:0<|a|<A^ — 1. We have thus proved the following 



Theorem 3.1 Suppose that rj G -"^(R") satisfies 

|r/(x)| < .4(1 + Ixl)--^, xGM", K>N + n, 

/n-|™.<». o<h<a'-i. 



and J-'q{Q) 7^ 0. Then the function 

N-l 

(/3- a) !(27ri) 1/5-^1 



|/3|=Oa</3 

satisfies the conditions (I14p . 

Suppose that r/(x) is radial, that is ??(x) = ip{r),r = \x\. Then d°'J^rj{0) = for any 
a = (ai, . . . , a„) containing at least one odd and we obtain the formula 



N-l 



.1) a27(^r?)-i(o) 



|/3|=o 27</3 V /; V ; 

Let us consider the special case r]{x) = vr'^/^e^l^l^ with J-ri{X) = e^'^^l'^'^. Denoting by 
Hp the Hermite polynomial of n variables defined by 

Hp{t)=e\'\\-dfe-\'\" 

we derive d'^ {J^rf)'^ {Q) = (— 7ri)l'''l//^(0). Note that -^^(O) = if 7 has odd components, 
otherwise H2^{^) = (-1)'^ (27) 1/7!. Hence 

«(-) = E '^M-)e-W' E (») 

|/3|=0 27</3 '• 

Assuming = 5|^|o and = 2M we find out 

M-l 



i=o |/3|=i 
= ^-n/2 ^ l_lLA^eH^'l' =^-"/24y_2)(|x|2)e-l-l' 



i=o 

where L^''^'* are the generalized Laguerre polynomials, which are defined by 

In this case we obtain the classical generating function ry(x) = L^(_^](|xp)e~l^l^ (see [9,13]). 

4 Examples 

Example 4.1 If N = 2, then formula (j20]) ^'zt'es 

7^(x) = ^-"/2/i _2 ^ (21) 



Figure 1: The graphs of {Ma — I) cosx when a = (sohd hne), a = 1/8 (dotted line) and 
a = 1/4 (dashed line). 



For the one- dimensional case the corresponding quasi- interpolant is 



X — km , _l^_h^\2 llh'^n 



Mau{x) = (ttP)"^/^ ^(^(fom) -/tVPan^(/tm))(l-2a )e 

For u{x) = cosx, x € M, the difference between u{x) and Mau{x) is plotted in{l\ by taking 
h = 0.1, T> = 2 for different values of a. 

Example 4.2 Now we are looking for quasi-interpolants of order 0{h'^). In the one- 
dimensional case 

= -ai, ^02^^ = fl? ~ "2, •^03^'* = + 2aia2 - as 
and formula [2D\) gives 

+ (4^-^) - l)x' + (8^3^) - 2A^-,''>)x% 

If ai = 02 = as = 0, then we get the classical generating function ri{x) = 7r^"'^/^e~^'^(3/2— 
x^). Assuming as = a2 = 0, ai = ±1/2 we obtain 

and the quasi-interpolant of order 0{h'^): 



MMx) = (ttV)-'/^ yiu{hm) ± ^v!{hm)){l ± ^^J^) ^-i-h^? /{h'v) 



With the choice oi = as = 0, a2 = —1/4 we obtain 7i{x) = vr ^/"^e and the quasi- 
interpolant of order 0{h'^): 

M2u{x) = (7rP)-i/2 Y^^^(^hm) - ^u" {hm)) e'^--^^^' '^^''^^ . 

mail 

In the Figures \^\^\^we show the error graphs for the approximation of u = cos(x) 
with the quasi-interpolants Mu, Miu, M2U, respectively. We have taken D = 2, h = 0.05, 
h = 0.1 and h = 0.2. The case of smaller h gives different pictures: it is clearly visible 
that the error oscillates very fast with Mu and Miu. 



If n > 1, by taking 02/3 = a and = otherwise, we obtain the radial generating 
function 

Haix) =(1 - ^ {a2f, + l/A){4x^f - 2))e-l-'lV-/2 
=(1 + n(2a + 1/2) - (1 + 4a)|x|2)e-l^lV-"/2 



for the approximate approximation of order 0{h'^ 

MaU{x 



(7rP)-"/2 {u{hm) + V a Au{hm))'Ha 



Example 4.3 Let n > 1 and set in ^ N = 2M and 

a-y = , if 7 has odd components, 



/X — km-. 



-)/!4l7l 



< hi < M - 1. 



Keeping in mind that 



^d'^''u{x) =A'uix) 

|7|=s 



we obtain 



M-l 



Q{{-hVV)d)u= 

j=0 



M-l 



4i 



h\=j i=o 

The function 7{{x) = 7r~"/^e~l^l^ satisfies conditions [T^ . In fact 

0, if a has odd components. 



n/2 



X e 



(-vri " 



(27)! 



i/ a = 27 . 



l7!22|7l' 

T/ien i/ie equations are valid for all P : < \P\ < M — 1 because of 



E«2(/3-a) _ 



a!22|°l (/?-a)!4l/3-"l a!4H ~ 

Therefore, a general approximation of order N = 2M is given by 

M-l 



0, ...,M-1. 



meZ" s=0 



2s 



s!4« 



-\x~hm\^/{h^V) 



If M = 2 then we find the 'fourth order formula" 



U(X} 



m) 



Au{hm)] e 



-\x-hm\'^/{h'^V) 



For M = 3 we obtain the "sixth order formula" 



(22) 



(23) 



(24) 



M(6)n(x) = (^P)-"/2 ^[n(/im) - ^Au{hm) + '^AMhm)]e~\''~^^\'/^'''^\ 



Note the additive structure of the formula 



\x-hm\^/{h^V) 



Example 4.4 We consider the second order partial differential operator 

n 

Bu= y^ bikdidku , 

i,k=l 

where the matrix B = {hik} £ M"^" is symmetric and positive definite. Define the quasi- 
interpolant 

Mu{x) = ^"^""^^TT y y B%{hm)e-^^-'^^-^"^^^^-^^^/^'^'^\ (25) 

Assume C such that B^^ = C^C. If we consider the linear transformation ^ = Cx and 
introduce f7(C) = u{x), in the new coordinates we have ( [12, p. 42]) 

Bu{x) = AU{^). 

Then ()25p will take the form 

^ ' meZ" s=0 

Keeping in mind (p3|) and ^ we rewrite 

MU{0= {-hVv)\('\d''U{0£/3{l,V) + Rh{^) 



=0 



where 



with 



and 



{hWy^'detC 
(7rP)"/2 



Now we use Poisson's summation formula on affine grids (see [11, p. 23]) 



yn/2 2^ 



' i 

■.2TT 



■ m (26) 



where we denote hy C = {C'^) ^ ■ In our case ri{x) = vr "/^e '^'^ and J^r]{\) = e '^^I'^l^ 
Since d^rr]{\) = {--kY Hs{-k \)e-''''\^\'' we obtain 



, if 6 has odd components, 

2. (-1)1"! (27)! 



vr 



7 



if6 = 2j. 



Formula (|26p applied to Sq, ^rii^es 



— j ^ a"J^?7(\/pC"'^i/)e2'^*<«'<^ a has odd components, 



2l°l7 V27r/ ^ ^ 



We deduce that 



E 

27</3 



(_l)l7l 



i \l/3-27l 



7!(/?-27)!4lTl V27r 



and 



vr 



I27-/9I 



7!(/3-27)!2l/5| , 

repeating the same arguments used in the proof of Theorem \2.1\ we derive that 
^|/3|ol < £1 for prescribed ei > and sufficiently large V, and the remainder is bounded by 

\Rh{0\ < CBihVvf'' P^Uho^ 

|/3|=2M 

with the constant cb independent ofU,h,T>. Hence 
\Mu{x) - u{x)\ = \MU{C) - U{C)\ 

2M-1 



=0 



=2M 



2M-1 



=2M 



where Cb depends only on the matrix B. Therefore the quasi-interpolant (j25p approximates 
u with the order 0{{h^/V)'^^) up to the saturation error. 



5 An application of formula (1241) 

Here we consider the approximation of harmonic functions. Suppose that Au = in some 
domain Q. C R". Then for any N = 2M and x G the Hermite quasi-interpolant (I24h 
has the simple form 

M(^)n(x) = Mu{x) = (vrP)-'"/2 ^ ^^^^^ ^-\,,-hm\^ / (^K^v) ^ ^^l) 



i.e., it coincides with the well known quasi-interpolation formula of second order. However, 
Theorem 12 . 1 1 indicates higher approximation rates. This will be studied here in more detail. 
First we consider the case = W\ Then 

and the series converges absolutely in M". Moreover, u has the analytic extension 

^(0 = f; ^ (C - x)^ CeC", (28) 

|/3|=0 

cf. e.g. [1,14]. Using formula (jlOp for the quasi-interpolant with the generating function 
we obtain 



Mu{x) - u{x) = £h,2M{x) + Tlh,2M{x) , (29) 

where 

meZ" 

2A/-1 , 
|/3|=1 meZ" a<f3 hVT> 

constitutes the saturation error and the remainder term has the form 

X E E C//3(-,Me-l^-'^-l'/(^'^) E "-^(^-^r (30) 

|/3|=2MmGZ" 0<a</3 



From ()22p we see that 



^ a! ^ (/3-27)! ^ 7! (/? - 27)! 22H 



which by using the representation of Hermite polynomials 



shows that 



Hence we obtain 

2Af-l 



£:.,2m(x) = >^ (- 



with the functions 



aMx,P) = (.Pr"/^ 5: ^/3(^)e-l^-"^l'/^, |/3| = 1,...,2M-1. 



(32) 



It follows from the definition of Hfi and from Poisson's summation formula that 

meZ"\{0} 

Thus the saturation error can be expressed as 

«M«w= E ^(-) 

|/3|=0 
mGZ"\{0} 

where 

|/3|=0 

is the Taylor polynomial of the analytic extension u. Note that p2p is valid for any M . 
Theorem 5.1 Suppose that the harmonic in M" function u is such that the series 

1/31=0 



converges absolutely for any y £ M". If yDh < 1, i/ien i/ze quasi-interpolant ()27p approx- 
imates u with 

Mu{x) - u{x) = lim £h,2M{x)= V u(x + 7ri/iPm)e-"'^l'^l%2-iK^>M, 

meZ"\{0} 

where u is the analytic extension of u onto C". 

Proof. We have to show that |7^/i,2Af (3^)1 ^ as M ^ 00. To estimate ()30p we rewrite 

^-1)^' 1 

1! (k - 

0<2j<k ■ 



r 

1)1 j Hk-iit)dt, 



2k 

which implies for \(3\ = 2M 



f3j>0 > i 



Consequently, the remainder TZh,2M takes the form 



(7rP)-"/2 ^ Up{x,hm) 

|/3|=2MmeK" 

z 

\x-hm\'^/{h'^V) 



xe 



where we use the notation Zj = {xj — hrrij) / {hyD). Then Cramer's inequahty for Hermite 
polynomials 

\Hk{x)\ < 2^/2^e^''/2 (34) 
(see [2,15]), leads to the estimate 

9(/3,+l)/2 / 

|/3|=2A/meK" l3j>0 Va/^j-I)! ^ 

< (^)'^(vrP)-"/' E E l^/3(a;,/im)|5;3(x-/im). (35) 

|/3|=2A/ meM" 

For the last inequality we use the notations 

/3j=0 /3,>0 ^V^P 



X 5; 5; |C/,(.,/.m)|e-l-MV(.^^) JJ / e*V2,, 



and the estimate 

e-'" J e''/^dt<\zj\e-'V^. 



By Q we obtain for harmonic u 

1 



oo 



|9'3+°w(, 

|a|=0 

which shows that 

|C/g(x, /im)| 5/3(x — /im) < 



E(FTR! E S,(x-/.m)|(x-/.m)"|. 



To get an upper bound of the last sum we write 

Sf3{x — m)\{x — m)" I 



m.6K" 



and note that for T> > Dq 



/3j>0mjeIR 



Hence we obtain the estimate 



5'/3(j; — m)|(j; — m)"| < 



cP(l"l+")/2 



2" J-J- V 2 



a,- + 2 



/3j>0 



with a constant c independent of a, (5, and P. Now we use that for Uj > 2 



< 



2"j/2 



which leads to 



/3i=0 /3j>0 



with another constant ci. Thus we derive from ([35 
I'^/i,2m(2;)I 



|/3|=2Af |a|=0 

oo 



|a|!|/3|!p(l°l+")/2 |a^+"u(x)| 
|/? + a|! a! 

|a|!|/3|!p(l"l+")/2 |a^+"n(x)| 
|/? + a|! 



n 2«./2 n 



■n 2-"^/2n J2(a,+1)A 





+ 1)! 




-1)! 



|/3|=2M|a|=0 ' V ^^^Q 

with a constant C not depending on M, V, h, and u. Since 



a|!|/3|! I{a + P)\ ^\a\mi^ + Pl^^ 



\l3 + a\l y a!/3! " \(3 + a\l a!/3! 
the remainder can be estimated by 



\n,,2M{^)\<c{hVv)'^ E E /r^75TT 

|^|=2Af|a|=0 V (« + /?)! 

Thus, |'7^/i,2Af(2;)| ^ for any fixed V if ([33]) holds. 



p{|"l+n)/2 -Q ^2(a,- + l)f5, . 

/3j>0 



Figure 5: The graphs of Mu — u with u{xi,X2) = e^^ cos X2, V = 2, h = 2 ^ (on the left) 
and h = 2"'' (on the right). 




Figure 6: The graphs of Mu — u with u{xi^X2) = e^^ cosx2, P = 3, /i = 2 ^ (on the left) 
and h = 2~'^ (on the right). 



Remark 5.1 The assertion of Theorem 15. H is a concrete realization of a general approx- 
imation result for analytic functions. Let u be an entire function in C" of order less than 
2. Theorem 7.1 in [7] states that the semi-discrete convolution 

uu{x)= Y: n^e-l^-^-l^/('^^^) 

with coefficients 

Um ■■= j e-^'^l^l' u{hm + mVhy) dy (36) 

differs from u by 

Uh{x) - u{x) = ^(^ + ^vr/iPm)e-'^'^l™l'e2'^*<""'^>/'' , 

meZ"\{0} 

(cf. also [8, Lemma 2.1]). It can be easily seen from ()28p and (j36p that the coefficients 
Um = (7rl?)~"/^u(/im) if the restriction of u to M"' is harmonic. 

We have applied the simple quasi-interpolant (j27|) to the harmonic function u(xi,X2) = 
e^i cos X2 in E? by assuming V = 2 (see Figure \5^, V = 3 (Figure [6|) , P = 4 (Figure [7|) , 
h = 2~^ and 2"''. The experiments confirm that the quasi-interpolation error Mu — u has 
reached its saturation bound also for large h because it does not decrease if h becomes 
smaller. 



Figure 7: The graphs of Mu — u with u{xi,X2) = e^^ cosx2, T> = A, h = 2 ^ (on the left) 
and h = 2~'^ (on the right). 



Let now u be harmonic in some convex domain Q. C and we consider the approxi- 
mant 



Mu{x) = {T:Vy'^ u{hm)e 



\x-hm\^/{h'^V) 



(37) 



Theorem 5.2 Suppose that the function u is harmonic in a convex domain 0, C M" and 
satisfies for a given N = 2M 



cu= Yl W^^^h^m n 



\f3\=2M 



Then for any e > and subdomain 17' C $7 there exists T> > and h > such that the 
quasi-interpolant (j37p provides for all x £ 0,' the estimate 



N-l 



U[X 



Mu{x)\ < C{hVvfcu + eYl {hVv)\l^\\d'^uix)\ , 



(38) 



where the constant C depends only on the space dimension. 
Proof. Analogously to the case = we obtain 

\u{x) - Mu{x)\ < \'Tlh,2M{x)\ + \£h,2M{x)\ 

with 

\£h,2M{x)\ < \u{x)\ (7rP)-"/2 J2 e-l^-^'"l'/('^'^) - 1 



2M-1 



+ (vrP)-"/^ Y 



hVv\m\d^u{x)\ 



X — km 



-\x-hm\^ /{h^V) 



and 



\nh,2M{x)\ < 



hVv\N 



|/3|=2M hmen 



see ([35]) . Since 

\Uf3{x, hm)\ < Wdf'uWL^^n) and ^ Sp{x - hm) < cP"/^ 
with a constant c depending only on n, we get the inequahty 

\nh,2M{x)\ < c{hVvfc^. 

To estimate \£h,2Mix)\ we use the functions cr^g given by (f3T]l and write 



Furthermore, for x G we derive 



< 6p{h-^ dist{x,dn),V)^ 



where 5^(r, P) , r > 0, denotes the rapidly decaying function 



6,{r,V) = sup {nV)-/' 



-\x-m\^/V 



|a;— m|>r 



Thus, for any domain 0, fixed parameter h and multiindex /3 we can find a subdomain 
^'l^ h ^ such that 

sup 5pih'Uistix,dn),V) < \\api-,V)\\L^ , (39) 

which gives 

|^.,,m(x)|<2''£'(^)'''^^K(-,P)||.. for all x G ' fl' O',,, . 

i/3i=o m=o 

Now we have to choose h such that W C which is possible since ^'f^ ^ as 

Remark 5.2 Obviously the assertion of Theorems \5.1\ and 15.^ can be extended to the 
case that a solution u of the second order equation 



hikdidku = 



i,fc=l 

in is approximated by the quasi-interpolant 

Mu{x)=^^^^-^ n(/im)e-<^-'(--'^™)'--'^-)/('^'^) 

with the matrix B = {hik}. 



6 Approximation of derivatives 

Here we study the approximation of derivatives using Hermite quasi-interpolation operator 
([5]). We introduce the continuous convolution (see [11]) 

Csv{x) = n Q{-6d)v{y)dy (40) 

where Q{t) is the polynomial in ([6]). 

Theorem 6.1 Suppose that 7i satisfies the decay condition ()13p with K > L + n, L G N, 
L > N. For any e > there exists V > such that for any function u G VF^(M"') 

L-l 

\Mu{x)-C,^uix)\<e Y,ihVvp\dMx)\+cAhVv)^ ^ P^^IIl^ , (41) 

|7|=0 |7|=L 



where the constant ci does not depend on u, h and \JT>. 

Proof. Suppose that the function u G W^(M"). The Taylor expansion ([8]) with N 
replaced by L gives the following form of the quasi-interpolant Mu in ([7]) 

\-n\ 



Mu(x) =Y.a^ ^^^'"^''' ^"+^u(x)a^(|,P,?^) 



|7|=0 |a|=0 



|t,|=0 |a|=L-|7| meZ" \ IIM U J \ tlM U / 



Similarly the Taylor expansion of u around y leads to 

N-\ i-l-|7l 



Ci,u{x) = V V (-(5)l"+'^l9"+^n(x)^ / T''n{T)dT 

|7|=0 |q|=0 

+ (-(5)^Va^ V ^/ r"W(T)i7a+^(x,x-r5)dT. 

|7|=0 |a|=L-|7| 

Setting 6 = /ia/P, we obtain 

Muix)-C^^u{x)=Y,c^j E ^^^^T d''+M^)£a{l,V,n) 

|7|=0 |a|=0 

1^1=0 |a|=L-|7| mGZ" ^ '^^''^ ^ 

- J r"W(r)[/a+-y(x,x - r(5)dr]. 



(42) 



hVV J 



Then 



|7|=0 |a|=0 

+ (/iVP)^^'K| 5] [||p^(,p, 7^)11^^ + / |r"H(r)|dT). 

I in I I r I I V Jr" J 

h\=0 |a|=L-|7| 

Proceeding as in the proof of Theorem 12.11 we deduce (|4ip . 



Theorem 6.2 IfTC satisfies the conditions ()13p with K > N + n and ()14p . then for any 

\Csv{x) - v{x)\ < 025"" Yl ll^^^ll^oo- (43) 

\a\=N 



Proof. The representation (|42p with L = N gives 



where the remainder .R^^-^ satisfies 

|i?^,^(x)|<5^ J] ^ / |r^W(r)|dr<C2<5^ (45) 

|Q|=Af J<a ' \a\=N 

If condition (fH|) holds, in view of (01]), we obtain (03]). 
Theorem 16.21 leads immediately to the next corollary. 

Corollary 6.1 Suppose that 7i satisfies conditions ()13p with K > N + n and (jl4p . T/ien 
/or suc/i that d^u^ W^{W), 

\C5d^u{x)-d^u{x)\<C25^ 1 1 9"+^^* 1 1 Loo- (46) 

\a\=N 

If the derivative d^Ti exists and satisfies the decay condition (jl3p then the convolution 
satisfies 

u{x) = (hVv) l^'C^g/j-^ ^(x) = Csd^u{x) (47) 

where 

Cs^QPn^x) = j^^ d^n Q{-Sd)u{y)dy 
and the quasi-interpolant Mu in ([7|) satisfies the equation 

d^Mu = {hVv)-\^\MQpj^u (48) 

with 

Mq,-^u = f^^pL\ Q (^-hVVdj u{hm). 

mgZ" \ V / 

Hence keeping in mind (j47p and (I48p we write the difference 

(/i\/p)l^l[9^Mu(x) - a'^u(x)] 

= (/i\/p)l'3|[a^Mn(2:) - d^Csu{x)] + (/i\/p)l^l [5^C5n(x) - d^u{x)] 
= [Mqpt^u{x) - Cs^Qp^u{x)] + {hVv)\'^\Csdf^u{x) - a^n(x)]. 

From (j46p and ()4ip we obtain the following result. 

Theorem 6.3 Suppose that 7i satisfies the conditions ()13p mi/i K > N + n and ()14p . 
Moreover suppose that the derivative d^H exists and satisfies the conditions (fT3]) wit/i 



K > L + n. The for any e > there exits V > such that, for any u G W^{W^) with 
L>N + \f3\, 

L-l 

|7|=0 

+c^{hVv)^-\^\ \\d'^u\\L^+C2ihVv)'' Wd'^^^uU^. 

|7|=L \a\=N 

We deduce that formula ([5]) gives the simultaneous approximation of the derivatives 
d^u with the saturation term e/i~I^L 
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